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ABSTRACT 

The reduction of structural vibrations in ships and submarines is a long standing 
concern of the Na\y. Waveguide absorbers are very elTcctive devices which can be ap- 
plied to this problem. This study evaluates the increase in vibration damping of a plate 
structure across a broad frequency range using light weight beam waveguide absorbers. 
Viscoelastic and constrained layer beam waveguide absorbers were studied both theore- 
tically and e.xperinientally. Impedances of the waveguide absorbers at the attachment 
point were predicted using both Bernoulli-Euler and Timoshenko beam theory for the 
viscoelastic beam and using sixth order beam theory for the constrained layer beam. The 
theoretically predicted impedances were compared with the experimental measurements. 
Results from random vibration tests of a plate structure showed significant increases in 
damping over a broad frequency range (100 Hz - 2 KHz) when the waveguide absorbers 
were attached on the plate. 
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I. INTRODUCTION 



Suppressing noise and vibrations of ships and submarines is very important for the 
Navy and has been one of the Naty's long-standing concerns. .Much work has been done 
on the development of vibration control means such as isolation, detuning, viscoelastic 
damping and dynamic absorption. All of these approaches have been studied rather 
thoroughly including their range of applicability and definite limits on their vibration 
and noise reduction capability within acceptable weight and volume increases. 

Recently, a waveguide absorber concept has been developed [Ref IJ. which can 
provide simple and highly effective vibration control over a wide frequency range. A 
"Waveguide" is a structure along which vibrational waves can travel. If one end of a 
waveguide is attached to a vibrating structure, some vibration energy will travel along 
the waveguide in wave forms. If the waveguide is treated with a high energy dissipation 
scheme, the damping of the waveguide causes the amplitude of the waves to decrease 
as they travel and the waveguide may be expected to remove vibrational energy from the 
structure. 

In previous studies by Ungar and Kurzweil [Ref 2] and by Ungar and Williams (Ref 
?], semi-infuiite beams and exponentially ta]'>ered semi-infmite beam waveguide absor- 
bers were studied. However, theoretical prediction of the driving point impedances 
showed wide discrepancy from the experimental results since the experiments were per- 
formed with finite length beams and the effects of viscoelastic materials on beams were 
not considered in theoretical calculations. 

In this investigation, two kinds of high damping beam waveguide absorber are 
subjected to study, viscoelastic beam and constrained layer beam waveguide absorbers. 
The driving point impedances are theoretically predicted for the finite length beams of 
highly viscoelastic behavior. Experimental studies are also performed on the impedances 
of these waveguide specimens and on their contribution to damping increase of a test 
plate structure. 
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II. THEORETICAL ANALYSIS 



A. LOSS FACTOR CONTRIBUTION OF A WAVEGUIDE ABSORBER 
Consider a linear structure. S, under a harmonic load 



F=fo^[noj[ ( 2 . 1 ) 

at location c (F igure 1 on page 3). The structure 5 will vibrate with the energy of vi- 
bration 






( 2 . 2 ) 



where M represents the total mass of the structure and is the magnitude of average 
velocity. The energy dissipated in 5 per cycle. depends on the original loss factor. 
of the structure and the following relation holds; 



For a lightly damped structure. is much smaller than 1 and is much smaller than 
27:11 Q. If a point a of a second body B. a waveguide absorber, is attached to body S 
at location b of body S . internal force. F*. will interact between body S and body B and 
some vibration energy. D will be dissipated in body D (Figure 1 on page 3). Moreover, 
this internal force inlluences the vibrational motion of body S and may reduce the en- 
ergy of vibration. IV. of body S under the same harmonic load at the same location c . 

1 he interaction between bodies 5 and D depends on the impedances of bodies .S and 
D at the attachment point, Z, and Zy Impedance is the ratio of the harmonic force. F. 
acting on a point a of the structure to the velocity. V, of the point. 

Z = -^. ( 2 .- 1 ) 

and can he describd Iw a complex number. 

Z=R + iX {2.5) 
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I 




Figure I. Vibrating structure uith waveguide absorber. 



where R is the real part of' Z and resistance component, or simply resistance, and X is 
the imaginary part of Z and reactive component, or reactance. 

The energy loss per cycle. D, of a structure by an attached waveguide absorber is, 
as shown by the previous study [Ref 2 . and Ref 3], 



D = 



^Rh 

0 ) 





( 2 . 6 ) 



where co denotes the vibration frequency in radian sec, represents the velocity of the 
structure at the attachment point of the structure before the attachment of the wave- 
guide absorbers. Zj and Z, denote the impedance of the waveguide absorber and the 
structure at the attachment point, and R^, represents the real part of the impedance Zy 
From the equations 2.2, 2.3 and 2.6, the energy loss factor. due to the 
contribution of an attached waveguide absorber can be expressed as follows; 
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Rh 



}j = 



coM 



11 + 






(2.7) 



Equation (2.6) shows that the vibration energy loss of a structure through the wav- 
eguide absorber depends on the magnitude of the velocity. 1', at the attachment point 
on the structure. If this attachment point is a node for which E) = 0. then the energy loss 
factor. )/. equals zero. If the attachment point is an antinode, then | -p- |- can be rela- 
tively large. Equation (2.7) indicates that the loss factor contribution is small if the im- 
pedance of the beam. , is considerably greater than the impedance of the structure. 
Z,. If Z* is much smaller than Z,. then the structure will not be alfected by the absorber 
and the vibration energy of the structure, II q in equation (2.2). will not be reduced much. 
Therefore, impedance matching between the structure and the waveguide absorber is 
important to get the highest energy loss factor, )/. So. analytic methods which can predict 
the impedance of beams using the theory of elastic wave propagation in beams are de- 
veloped in the next section. 

B. IMPEDANCE OF BEAM WAVEGUIDE ABSORBERS FROM WAVE 
PROPAGATION THEORY 

1. VISCOELASTIC BEAM USING BERNOULLI-EULER BEAM THEORY 
Previous research [Ref 3] considers the impedances of the infinite viscoelastic 
beam using Bernoulli-Euler beam theoiy. In this section, this previous work is expanded 
to include the finite viscoelastic beam. The viscoelastic beam will be excited harmonically 
at its center by a transverse force, F as shown in Figure 2 on page 6. For a beam which 
is free of lateral loading and under an assumption that cross-sectional areas remain plane 
and normal to the neutral axis, the equation of motion bccames [Ref 2). 









~ + pA r = 0 

c.v dr 



( 2 . 8 ) 



where y(.v.f) is the transverse displacement, and E’ is the complex moduli of the beam: 

r = £U + /;/J. (2.9) 



and. 

E : modulus of elasticity of the beam 
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: loss factor of the viscoelastic material 
m : mass per unit length of the beam 
p : density of the beam material 
A : cross sectional area of the beam 
/ : moment of inertia of the cross sectional area 

'Ihe general form ofy(.Y.r) for a propagating harmonic wave in this Bernoulli- 
Euler beam is 






r M'x-a,t) 
I 



(2.1U) 



where A' is the complex wave number, and o) is the frequency of the propagating wave. 
Inserting the wave equation (2.10) into the equation of motion (equation 2.8). the dis- 
persion relation is obtained as follow: 



£V(A"V = pA(x? 



( 2 . 11 ) 



This equation gives the four dilferent complex wave numbers as function of frequencies. 

V / ~ _iM tan~'«,,„ tan"’;/.,. 

1 + ’he ) %os{ ) - i sm( )j (2. 1 2.nj 



/<2 = ~‘K\i 1 + ’lie ) ^'\cos{ ) - i sin( "" )] (2. 1 2.b) 



tan ? 7 ,, 



, ^’h-e 



, -> , -> 

V / 7^ -\id tan tan n..„ 

f<3 =- /<(v ’I've) ' [cns( ^ ) - / sm( )] 



(2.12.C) 



_ 1 1 
k: = kQ 1 -b )-''‘^(cn5( ^ ) - / sin( )] 



(2.12.cf) 



where subscript ve is an abbreviation of Viscoelastic beam using Bernoulli-Euler beam 
theory and k is the wave number, given by 



' El El ' ’-Ci 



(2.13) 



) and ty denotes the 



/ T 

r denotes the radius of gyration of the cross sectional area, {\J ~ 
longitudinal wavespeed in the beam material (\~^)- Therefore, the equation (2.10) is 
expanded using wave number /<h kl, ki and k^. 
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Y\e“'^y 



e “ 



i!oX\ — iVof 



(2.14) 



where F,, and Fj are constant.s whose values depend on the boundary conditions 
and X is the axial coordinate. The complex transverse force, F, is denoted by 



F(.r) = 



d\iU) 

cix^ 



(2.15) 



a. INFINITE VISCOELASTIC BEAM 

Let's consider tlie infinite viscoelastic beam. If rf.r) is to remain finite for a 
large .v, then F^ must vanish. If there is no wave coming toward the origin from the po- 
sitive X direction, then must vanish. F, and in equation (2.14) are determined from 
two boundary conditions at .r = 0, as functions of y and B, as follows: 

y(0)=y (--16) 
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(2.17) 



dx 



= 6 



At the driving point, the input transverse Ibrcc. r(0). can be calculated using 
the equation (2.15), which can be expressed using impedance and velocit}' at .v = 0 as 
follows: 



f(0) = Z, F(0) + Z2Q(0) (2. 1 8) 

The complex transverse velocity, V, and the complex angular velocity, 12. are related to 
the corresponding displacement, y and 6 of equations (2.16 and 2.17). as follows: 

^ = = (2.19) 

Z, represents the impedances of the transverse force at the driving point ol' the center 
of the infinite viscoelastic beam. If there is no rotation at the center of the beam, that 
is Q(0) = 0. the equation (2.18) becomes 

T(0) = Z,r(0) (2.2U) 

Therefore, the impedance of an infinite Bernoulli-Euler beam at the dri\ing point. Z,„,, 
can be calculated from 



ZireiO) = 



^vg(b) 1 
ioj J' 



( 2 . 21 ) 



where subscript ive is an abbreviation of the Inlinite Viscoelastic beam using Beinoul- 
li-Euler beam theory. 

b. FINITE VISCOELASTIC BEAM 

For a finite beam, four constants T,, \\ and F,, are determined iVom four 

boundary conditions at .v = 0 and .v = /, as follows: 

T(0)=y (2.22) 



dx 



= 0 



-E'^I 



d\v{[) 

dx^ 



= 0 



(2.24) 
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(2.:5i 



E^l 



d\xil) 



= 0 



From the same considerations as those in the previous section, the impe- 
dance 



ZUO) = 



io) 



(2.26) 



is determined as function of frequencies, where the subscript fve is the abbreviation of 
the Finite Viscoelastic beam using Bernoulli-Euler beam tlieoiw. 

2. FINITE VISCOELASTIC BEAM USING TIMOSHENKO BEAM THEORY 
The Bernoulli-Euler beam theor\' was obtained using two assumptions which 
neglect shear deformation and rotational inertia. The Bernoulli-Euler beam theory may 
be useful for low frequencies and long wavelengths. However, the Timoshenko beam 
theory gives more accurate simulation at wider frequency ranges and wave lengths. 
Therefore we will analyze the impedances of the finite viscoelastic beam using Timosh- 
enko beam theory in this section. Two independent variables, the transverse displace- 
ment, j(.v.r), and angular displacement, ^f.v./). are defined in Timoshenko beam theory. 
Beam motions are governed by the following equations (Ref. 4|. 






-2 

C 



cx 



C(p 

cx 



c r 



cr 



(2.27) 



2 C^d) , C^(f) 

C1-) "h A'q ( 0) 



i, 

cx 



cx 



(2.28) 



ct 



1 , kG^ 

where a\ — — a; 



— and A'^ = ^ The complex moduli G^ and E' are defined as: 






(2.29) 



E = E[ \ + ii} 



E30) 



where subscript vt is an abbreviation of Viscoelastic beam using Timoshenko beam the- 
ory, and 

G : shear modulus of the beam 

E : elastic modulus ol'the beam 

>!,, ; loss factor of the viscoelastic material 
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A : cross sectional area of the beam 
/ : moment of the inertia of the cross sectional area 
p ; density of the beam material 
K ; shear deformation factor 

The general forms yl.v.r) and 4>{x,i) for the propagating harmonic wave in a Ti- 
moshenko beam are 



y(.r,r) = (2.31) 

(2.32) 



where k'^ is the comple.x wave number, and co is the frequency of the propagating wave. 
Inserting the wave equations (2.31 and 2.32) into the governing beam equations (2.27 
and 2.28) becomes 



(o)^ - «r(A-")^) To + ( -/c7?A")Oo = 0 



(2.33) 



{ia^k^k^) To + (r^- - afAj - a2'(A^j-)cho = 0 



(2.34) 



The dispersion relation of waves is obtained from equations (2.33 and 2.34) 
a^aj{k^f — (af 4- ci2)oj~{k^)^ -b {o/ — a^k^co^) = 0 
with the relationship between ip and d>p. 



On = 



po) 






ixO^k 



>0 = /^^', 



where 



R,= 



pco" — kG’^(A^)^ 



(2.35) 



(2.36) 



(2.37) 



The equation (2.35) gives the four different wave numbers as function of wave frequen- 
cies. 



Ar = 



(U] -l- Ui T 



2s. 



y j I'll 

(c?f -t- Ujj'co — aCctiiiT — aCAoOJ-) 



2«r«* 



(2.38.0) 
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^2 =- 



k; = 



k:=- 



2 "^2 "'"'4 *'"^4 ''0'’ 

(67| + + \ H" <^^V'co - i7fc?;f(t) — czf/<QC'.r) 

-> 2 2 
2ajc?2 



(c7|“ + (-h)^^ ~ \ (^^i~ cxjy'co^ — a^al(co^ — a^k^]o)~] 



2 2 2 ' '^"'4 0 ''4 ') 1 "? 

(af + a2)co^ — Y ia\ + c?2)"tu — afai (ty — «!■/:,) cy') 

-,2 2 
2a, (72 



12.38./)) 



(2.38.C) 



(2.38.(y) 



Therefore, the equations (2.31 and 2.32) are expanded using Af, A|. A) and kl 

! .1 /I' 1 '^T’c , I- Ik^x 



The equation (2.40) becomes the following equation 



where 



Rn = 



pco^ — KG^ik^y 
iKG'^k^ 



for «= 1. 2, 3. 4. 

For a finite viscoelastic beam with a leng 
at .V = 0, four boundary conditions are defined as: 

>■(<•>) =y 

4>{0) = 0 



A/(/) = 

cx 



ov 



vx 





(2.39) 




(2.40) 




(2.41) 




(2.42) 


which is free at j 


(• = / and excited 




(2.43) 




(2.44) 


1 


(2.45) 


= 0 


(2.46) 



From those boundary conditions, the constants 1’,. 13 and Iq are determined 
as functions of y and <i). At the driving point, .v = 0, input transverse force, 
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/"(O) = — S(0), can be calculated using equation (2.46), which can be expressed using 
impedances and velocity at the driving point as follows: 

/••(0) = Z, (■(()) + Z2(b(U) (Z47) 



where the complex transverse velocity. 1", and the complex angular velocity. <l>, are re- 
lated to the corresponding displacement as follows: 



.t’ (j) 



— i(o 



(2.48) 



If there is no angular displacement, 4>, at the center of the beam, that is, <h(b) = 0. the 
equation (2.47) becomes 



F{0} = Z, f '(0) 



(2.49) 



Therefore, the impedance of the Timoshenko beam at the driving point, Z^,,, can be cal- 
culated using the following equation. 



Z/„(0) = 



frJO) 1 

—i(x> y 



(2.50) 



where subscript /vr is an abbreviation of the Finite Viscoelastic beam using Timoshenko 
beam theory. 

3. FINITE CONSTRAINED LAYER BEAM 

In this section, impedances of the constrained layer beam, which has viscoelastic 
material between the two elastic layers, are studied. The theory of a damped sandwich 
beam was developed by R..A. Ditaranto [Ref. 5], who extended Kerwin's [Ref 6 ] bas- 
ically similar analysis. D.J. Mead and S. Markus [Ref 7] expanded the concept by con- 
sidering these earlier works. The result of their elforts is a sixth order differential 
equation of motion which is expressed in terms of the transverse displacement, y, and 
longitudinal displacement, u. for the constrained layer beam. The constrained layer 
beam will be excited harmonically at its center by a shear force, S, as shown in 
Figure 3 on page 12. This equation of motion leads to the impedances of the sandwich 
beams with six boundry conditions, which are used for evaluating impedances depending 
on frequency, viscoelastic material properity, and two elastic material properties. 

From the [Ref 7], the governing sixth order dilTerential equation is 
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^6 ^4 ^2 

-rv-s»(‘+2)T^ = -^iTT-'P^i '-'51) 

CX CX ‘-^J OX 

where 



c i- 

^ = -m 

ct' 



(2.52) 



with the relationship between _v(.v.r) and u{x,t) expressed as: 



d'u r, c\v 

- B,,u = -rr- 



(2.53) 



The pctrameters in this equation are defined as follows: 






+ ■ 



H. '£,//, £j//3 



■) 



I2.54.cr) 
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^ C,^ E, 7 /,E 3//3 

Dj ^ £,//, + £32/3 ^ 


(2.54.6) 


E\I^ 4- £ 3/3 


(2.54.C) 


Q =//2 + ^(//, +7/3) 


(2.54..f) 




(2.54.C) 




h] 


(2.54/) 


'^=-rr 


G*' = G(l + inc) 


(2.5.'^) 



where subscript c is an abbreviation of constrained layer beam, and 
G ; shear modulus of the viscoelastic layer 
E : elastic modulus of the elastic layer 
: loss factor of the viscoelastic layer 
/ : moment of the inertia of the cross sectional area 
H : thickness of the constrained layer beam 
m : mass per unit length of the whole three layer section. 

The general form of transverse displacement, and longitudinal displace- 
ment, for the propagating harmonic waves in a constrained layer beam are 

(2.56) 

u{x,t) = (2.57) 



where A." is a complex wave number and w is the frequency of the propagating wave. 
After inserting the wave equations (2.56 and 2.57) into the sixth order beam equation 
(2.51), the dispersion relation is obtained as follows: 






tno) 

D, 






Dt^nuo" 

~dT 



= 0 



(2.58) 
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There are six different roots with six dilTerent eoinplex wave numbers in the above 
equation, which means a wave of six different modes propagate through the constrainei 
layer beam. From the equations (2.56 and 2.57), we may rearrange equation (2.53) 

Uo = RoYo r-59) 



wliere 






(2.60) 



The general form of the transverse displacement equation (2.56). y. and longi- 
tudinal displacement equation (2.57). u , are obtained using the six complex wave num- 
bers and a ratio Rf^. The equations (2.56, 2.57 and 2.59) as follow; 



/ ^ U<.y X , T- lk-,X , T' IKtX , TT 

r(,Vj) = (} |<^ ^ +^ 2 ^ ' + ^4^ 



ik^x 






(2.6i) 



a(.v,0 = (/?,f>'^ 



+ Rj + Ri + Ra ^ 



-b 



n r • ik<X , n T - — ;coJ 

7^5 fjC - ^ )e 



(2.62) 



where 



R 



n 



iRbZOX 

{{O- +D,)E,II,C, 



(2.63) 



n = 1,2, 3, 4, 5, 6 For a constrained layer beam with a finite length / which is free 
at .V = / and excited at one end at .v = 0, the six boundary conditions are defined as fol- 
lows [Ref 7j: 



_V(0) = y 


(2.64) 


- 1) 

ClX 


(2.65) 


»(0) = u 


(2.66) 



S(f) = A 



c y 
cx 



- 



cy 

dx 



, D,B,Z ^ £■,//,+ 7f3//3 . 

+ — — ( F77 )u = 0 



C 



EiH 



(2.67) 
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m = D, 



( 2 . 68 ) 



"" r-QE,ii, 4^ = 0 

• cx 



CX 



Ml) = EMC ^ 

C.Y 



(2.69) 



The six constants (/ =1, 2. 3, 4, 5 and 6) are determined from tliese six boun- 
dar>’ conditions and are a function of j, d and ii. At the driving point, .v = 0. the input 
transverse force. /'(()) = — 5(0). moment. 3/(0), and normal force. A’(0). can bo calculated 
using equations (2.67, 2.68 and 2.69). These input force and moment can be expressed 
using impedances and velocities and angular velocity at the driving point as (ollows: 

F(0) = Z, I '(0) + 22^(0) + Z3 L'(0) ( 2. 70) 

Where the complex transverse velocity, F, and the complex angular velocity, 
Q. and the complex longitudinal velocity. U, are related to the corresponding displace- 
ment as follows: 



- ^ ^ - 
y n 



(2.71) 



If there is no rotational displacement, 9, and no longitudinal displacement, u. at the 
center of the beam, Z[ represents llie impedances of the transverse force at the driving 
point of the constrained layer beam. That is, £2(0) = 0 and L'(0) = 0. Therefore, the 
equation (2.70) becomes 



F(0) = ZiF(0) (2.72) 

Therefore, the impedances of the constrained layer beam at the driving point, Z^, can 
be calculated using the following equation. 



ZciO) = 



m 1 

—i(o y 



(2.73) 



Where subscript c is an abbreviation of the constrained layer beam. 

C. COMPUTER IMPLEMENTATION 

The computer programs using Reduce and I'ORTRAN languages are used for cal- 
culation of the theoretical impedances from the beam theories in this paper. Reduce was 
developed by Anthony C. Hearn [Ref 8], and is a programming language which solves 
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algebraic operations non-nuinerically. It can manipulate polynomials in a varicy of forms 
(substitution, differentiation and integration) and solves one or more simultaneous al- 
gebraic equations. After the impedances were evaluated by a Reduce program with 
constant variables, a FORTR.'W program calculated the impedance with real dimen- 
sions and properties of the viscoelastic and constrained layer beam waveguide absorber 
for the 0-2000Hz frequency range. All constant variables of the Reduce program were 
delined in each FORTR.AN program. We tried to use the same notations for the Reduce 
program, the FORTR.AN program and the te.xt (Appendix A). 

Three Reduce programs and a FORTIT.AN program were used to study each theory. 
The first Reduce program delined the moment, force and displacements which included 
the constants F,,. where n equals 1 and 3 for an infinite viscoelastic beam; n equals 1 , 2, 
3 and ’ for a finite viscoelastic beam and n equals 1, 2, 3, 4, 5, and 6 for a constrained 
layer beam. F„ can be solved with boundary conditions. After running the first Reduce 
program, new constant terms are defined for The second Reduce program includes 
those new constant variables which were the results of the first Reduce program. The 
second Reduce program gives the constant variables Y„. Finally, the third R.educe pro- 
gram can evaluate the impedance of the beam with longitudinal and angular displace- 
ment (and normal displacement for the constrained layer beam). The process of the 
constrained layer beam is similar up to the second Reduce program. However, the 
computer finds it hard to handle 6 simultaneous equations with 6 non-numcrical vari- 
ables, Therefore, we solved Tj, ); and in the second Reduce program and then 
substitute these values to and in the third Reduce program. 

Wave numbers as functions of frequencies. Le. dispersion relation, of waves in wav- 
eguide absorbers were described in previous sections. It was easy to find the wave num- 
ber of a viscoelastic beam through Bernoulli-Fuler beam theorx- and Timoshenko beam 
theory in equations (2.12 and 2.38). However, wave numbers of the constrained layer 
beam were expressed with a complex six order polynomial. Thcrelbre, they were calcu- 
lated using Newton's method. The equation (2. 58) is rearranged using real constants 
h,, c, and ^4 ^'id imaginary constants h, and c(. 

+ (b, + ibi)[K'^f + c>(A'’^) -b {d, + u/i) = 0 (2.74) 



where 



K"" = ( 4 ’^)- 



(2.75.U) 
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b,= Re{R,+ B,Z) 


{2.1 5. b) 


b,= lm{B^ + B^Z) 


[2.15.C) 


—fUCO" 


(2. 75.(7) 


d,-Re{ ) 


(2.75.C) 




(2.75,/) 



Let the complex variable K’' = i + iq (i and q are real numbers). The equation {2.74) 
becomes two equations coupled with two variables (/ and q). 

J{t.q) = P + b,.P + ( -3(7“ - 2biq + + ( -l\q^ + 4) (2.76) 

g[t,q) = q^ + biq^ + ( -3p - 2b, t - c^q + ( -b/ - (/,-) {2.11) 

Starting from initial guesses (t,, ^o) which are near the roots iteration was continued until 
either the successive t and q values converged to the certain values or the value of the 
funetion is sulTiciently near zero. After a few iterations of equations (2.76 and 2.77) 
using a FORTRAN program (Appendix A), we can get three different square of wave 
number values. A'% were obtained. So, finally, this method caleulates six different wave 
numbers. from ± K' , where n equals 1. 2, 3, 4, 5, and 6. 

Impedances of a 20" viscoelastic beam were calculated using these eomputer pro- 
grams for viseoelastic beam prop)crties of 200.000 psi constant shear modulus and 0.2 
energy loss factor. The results of these simulations are shown in Figure 4 on page 19 
and Figure 5 on page 20 which reflect sharp resonance frequencies. These sharp reso- 
nance phenomena reduced for the same beam with higher energy loss factor of 0.5. For 
eases using the Bernoulli-Euler beam theory, the impedances of the finite ^ iscoelastie 
beam approached to the infinite viscoelastic beam waveguide absorber with increasing 
energy loss factor and with increasing frequency. The impedance using Timoshenko 
beam theory differs from the impedance from the Bernoulli-Euler beam theory 
(Figure 6 on page 21 and Figure 7 on page 22), because Timoshenko beam tJieory 
considers shear deformation and rotary inertia effect of the viscoelastic beam. 
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Impedances of a 20" constrained layer beam with viscoelastic layer properties of 200 psi 
constant shear modulus and 0.5 and 1.0 energy loss factors were calculated. Figure S 
on page 23 and Figure 9 on page 24 compare the impedances vs. frequency depending 
on the energy loss factor of the constrained layer beam. The impedance of the high 
damping constrained layer beam shows the smooth and shifted resonance frequencies 
compared with the low damping constrained layer beam. 
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Figure 4. Real part of the 16 " viscoelastic beam uaveguide absorber theoretical im- 
pedances (?/ = 0.2) . 
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Figure 5. Imaginary part of the 16" viscoelastic beam waveguide absorber theore- 
tical impedances (;/ = 0.2) . 
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Figure 6. Real part of the 16 " viscoelastic beam waveguide absorber theoretical im- 
pedances (?/ = 0.5) . 
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Figure 7. Imaginary part of the 16" viscoelastic beam waveguide absorber theore- 
tical impedances (?/ = 0.3) . 
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Figure 8. 



Real part of the 20" constrained layer beam uaveguide absorber theore- 
tical impedances. 
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Figure 9. 



Imaginary part of the 2U" constrained layer beam naveguide absorber 



theoretical impedances. 
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HI. EXPERIMENT 



A. IMPEDANCES OF THE TEST PLATE 

A rectangular aluminum plate with clamped boundary condition was selected as the 
test structure, because it has enough closely space modes in the frequency range 
( 100Hz-2Kz). An aluminum block (30" x 22" x 2") was carved inside to ntakc a 24" x 16" 
X 0.3125" plate with clamp boundar\' condition. 

The impedances of the test plate were measured from impact hammer tests. The 
aluminum plate was excited by a PCB 0S6B03 impact hanuner. ln|)ut impact force was 
measured by a force transducer at the tip of the impact hanuner. Responses were meas- 
ured by an Endevco 2250A-I0 piezoelectric accelerometer which was attached under the 
impact point for locations 1, 2 and 3 (Figure 10 on page 26 and Figure 11 on page 
27). Location 1 is the center of the plate. Location 2 is the symmetric point of location 
4, where a shaker will be attached during the damping measurement tests. Location 3 is 
a general point which is not on any line of symmetrex A baseband measurement for each 
location was taken from 0 to 2000Hz. An average of 15 data samples was taken for each 
measurement and analyzed using a HP-3562A dynamic signal analyzer, impedance of 
the aluminum plate at each location was determined after a predetermined data 
processing procedure using the HP-3562A dynamic signal analyzer. 



25 




Figure 10. Arrangement for test plate impedance measuieinent. 
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Figure 11. Test plate impedance measurement configuration. 





B. IMPEDANCES OF THE WAVEGUIDE ABSORBER 



Two viscoelastic beams and two constrained layer beams were selected as the wav- 
eguide absorber specimens. Dimensions and physical piojierties are shown in Table I. 
Mechanical properties of a viscoelastic material depend on temperature and vibration 
frequencies. The loss lactors and shear modulli of the viscoelastic material used in the 
specimens at room temperature {75'^F) are shown in Figure 12 on page 30 and 
Figure 13 on page 31. These data were obtained from Design Data Sheets in a book by 
N.ishif. et. al. [Ref. 9]. .Analytic expressions of the loss factors and shear moduli were 
determined from curve fits to these data. The following four equations represent the loss 
factors and shear moduli of the LD-400 and lSD-112 viscoelastic material. 



[-0.527321 



log(- 






(3.1) 



;/^= 1 . 1(3 



[-0.39061 log(-^) 



, 1 - 53 , 



(3.2) 



Gv = 0.00002503/ -0.1752/ +457.5S83/-I-29280 



= -0.00001258/ +0.2472/+ 74.988 



(3.3) 

(3.4) 



f represents the frequency, t/, and represent the loss factor of LD-400 and lSD-112 
viscoelastic material. G, and G, represent the shear modulus of LD-400 and ISD-112 
viscoelastic material. Figure 14 on page 32 shows waveguide absorbers which were used 
in this study. 

Impedances of waveguide absorbers at the attachment points were determined horn 
steady state random vibration tests. The waveguide absorber was excited by a Wilcoxon 
F4, F7 vibration generator using a band limited white noise signal (0- 2000Hz) from a 
HP-3562A dynamic signal analyzer (Figure 15 on page 33). The waveguide absorber was 
mounted on the Wilcoxon F4 F7 vibration generator (Figure 16 on page 34). Force and 
ao'eleration were measured using a force transducer and a piezoelectric accelerometer in 
the impedance head of the Wilcoxon F4 F7 shaker. These .signals were analyzed by a 
HP-3562A dynamic signal analyzer. 100 data samples were averaged for each measure- 
ment. 

The measurement of impedances included the impedances of bolts and nuts and 
connecting aluminum pieces. Therefore, the impedances of bolts, nuts, and connecting 
aluminum pieces, the so-called mass elTect, were measured separately. Total impedances 
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minus mass edect represented the impedances of the waveguide absorbers. Tiiese 
experimental results are compared with the theoretical predictions in Chapter IV. 



Table 1. DIMENSIONS AND PROPERTIES OF THE WAVEGUIDE ABSOR- 



BER SPECIMENS. 





SPECI- 
.MEN I 


SPECI- 
MEN 2 


SPECI- 
MEN 3 


SPECI- 
.MEN 4 


WAVEGUIDE ABSOR- 
BER TYPE 


VISCO- 

ELASTIC 


VISCO- 

ELASTIC 


CON- 

STPv.AINED 

LAYER 


CON- 

STRAINED 

LAYER 


MATE- 

RIAL 


VISCO- 

ELASTIC 


LD-400 


LD-400 


ISD-II2 


ISD-112 


CON- 
ST FUA I NED 
LAYER 






A1 


Al 


LENGTH {in) 


16.0 


20.0 


16.0 


20. 0 


WIDTH {in} 


1.0 


1.0 


1.0 


1.0 


THICK- 

NESS 

(m) 


VISCO- 

ELASTIC 


0.35 


0.35 


0.065 


0.065 


CON- 

STIUAINED 

LAYER 






0.0625 


0.0625 


DEN- 

SITY 

^ in^ ’ 


VISCO- 

ELASTIC 


0.0550.5 


0.05505 


0.03663 


0.03663 


CON- 

ST1UVINED 

LAYER 






0.09832 


0.09832 
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Figure 12. 



Loss factor vs. frequency for a LD-4U0 and a ISD-112. 
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Figure 13. Shear modulus vs. frequency for a LD-40U and a ISD-112. 
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Figure 15. Arrangement for waveguide absorber impedance measurement. 
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Figure 16. 



VVilcoxon F4/F7 shaker with ^^aveguide absorber mounted. 
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C. DAMPING OF THE PLATE 

Modal damping values of the plate were determined from structural transfer func- 
tions at a driving point and between input and output points on the test plate which 
were measured using random excitations. The MP-3562.\ dynamic signal analyzer was 
used in the measurement of the transfer functions. The Wilcoxon F3 vibration generator 
excited the plate at location 4 using random oscillation signals from a MP-3562.'\ dy- 
namic signal analyzer. Input forces and output responses at the driving point were 
measured using a force transducer and a piezoelectric accelerometer in the impedances 
head of Wilcoxon F3 shaker (Figure 17 on page 36). The damping measurement con- 
figuration is shown in Figure IS on page 37. It was necessary to measure responses at 
several points other than the driving point to accurately measure damping values of 
closely spaced modes. These were measured using the Endevco 2250A-10 piezoelectric 
accelerometer. .Modal damping values of all the modes in the frequency range of 
100- 2000Hz were determined from the response functions using the Half power Band 
Width method [Ref. 10] and a curvefit method [Ref. 1 1] 

Modal damping values of the plate without a waveguide absorber were measured 
as the baseline. A waveguide absorber was attached to the plate at location 1. The 
waveguide absorber was changed after each measurement of the damping for four dif- 
ferent waveguide absorbers which were described in section 4-B. After the damping was 
measured at location 1, measurements were taken at location 2 and location 3. 

The most effective waveguide absorber at each location was selected. The 20" vis- 
coelastic beam waveguide absorber, the 16" constrained layer beam waveguide absorber 
and the 16" viscoelastic beam waveguide absorber were most effective at location 1, 2 
and 3. respectively. Therefore, the damping with two waveguide absorbers attached at 
locations 1 and 2, and the damping with three waveguide absorbers at locations 1. 2 and 
3 were also measure, to see the effect of multiple waveguide absorbers. 
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Figure 17. Arrangement for test plate damping measurement. 
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Figure 18. Damping measurement configuration. 
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IV. RESULTS AND DISCUSSIONS 



A. IMPEDANCES OF THE TEST PLATE 

Impedances of the test plate were measured from input hammer tests as described 
in Section 3. A. Since the aluminum test plate is a lightly damped structure, the impe- 
dances at resonance frequencies are very low and at anti-resonance frequencies are very 
high. Phases with 180'-^ phase shift are near + 90^ Iroth at the resonance and at the 
anti-resonance frequencies. 

Since locution 1 is the center point of the plate, only small numbers of resonance 
and anti-resonance frequencies appear {Figure 19 on page 39 and Figure 20 oir page 
40). Neither location 2 nor location 3 are on any line of symmetry of the test plate; all 
the resonance and anti-resonance points appear in the impedance vs. frequency dia- 
grams. (Figure 21 on page 41. Figure 22 on page 42. Figure 23 on page 43 and 
Figure 24 on page 44). However, the anti-resonance frequencies at location 2 and those 
at location 3 are dilferent. 
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Figure 24. 



Impedance phase of test plate at location 3. 
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B. IMPEDANCES OF THE WAVEGUIDE ABSORBER 

Two difFerent kinds of waveguide absorbers, a viscoelastic beam and a constrained 
layer beam, are investigated in this study. In a previous study by L'ltgar and Williams 
[Ref. 3|. impedances were predicted only for the infinite viscoelastic beam using Ber- 
noulli-Euler beam theory and the theoretical prediction showed wide difference Ifom the 
experimental results. In this study, this previous work is extended to include the impe- 
dances for the finite viscoelastic beam. Since Bernoulli-Euler beam theory does not in- 
clude the shear deformation and rotary inertia effect of the viscoelastic beam, the 
predictions from the Timoshenko beam theory are studied to investigate these effects on 
wave propagation and impedances, and are compared with the results from the Bcr- 
noulli-Euler beam theory. The impedances of the constrained layer beam waveguide ab- 
sorber were predicted using the sixth order beam theory. 

The impedances from these theoretical predictions are compared with the impe- 
dances from the experimental measurements for the 16" and 20" length waveguide ab- 
sorbers of two different kinds. In the calculation of theoretical predictions, the loss 
factors and shear moduli of viscoelastic material were varied as function of frequency 
as described in Section 3.B. During the experiment, the Wilcoxon F4 F7 vibration gen- 
erator produced characteristic errors which occured at different frequencies depending 
on the mass of the waveguide absorbers. Therefore the viscoelastic beam has unreliable 
data around 700Hz and the constrained layer beam has unreliable data around 3.501 Iz. 

Figure 25 on page 47 and Figure 20 on page 4S show the three theoretical impe- 
dances compared with experimental data. I'he Bernoulli-Euler beam theory for finite 
beams shows some standing wave effect at low frequency range due to reflection at the 
free end of the beam, but approaches the results of the infinite beam as frequency in- 
creases. The results of the Bernoulli-Euler beam theory do not quite correspond to ex- 
perimental I'esults for the viscoelastic beam. Impedances from the fimoshenko beam 
theory closely follow the experimental results for viscoelastic beams as shown in 
Figure 25 on page 47 and Figure 26 on page 48. These show that the shear deformation 
and rotary inertia effects are significant for the impedance prediction of beam type 
waveguide absorbers. 

Figure 27 on page 49 and Figure 28 on page 50 represent the real and imaginary 
parts of the 16" constrained layer beam waveguide absorber and Figure 29 on page 51 
and Figure 30 on page 52 represent the real and imaginary parts of the 20" constrained 
layer beam waveguide absorber. In both beams, the theoretical i'esults of the impedances 
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ii'iing the sixth order beam theory and the resonance frequency patterns are almost the 
same as the experimental result. 

Both Timoshenko beam theory and the sixth order beam theory predict resonance 
frequencies shifted amount from the experimental resulrs as frequency increases, and the 
impedance magnitudes of experimental results are bigger than the theoretical results. 
These diderence between theoretical and experimental impedances may derive from the 
inaccuracy of shear modulus and loss factors which values were selected from 75'"F room 
temperature and changed depending on frequency. Mowever. the theoretical prediction 
using Timoshenko beam theory for the viscoelastic beam and the sixth order beam the- 
ory for the constrained layer beam follow the same trend (Figure 12 on page 30 and 
Figure 13 on page 31). 

Impedances of waveguide absorbers are very low and do not .show any sharp vari- 
ation for all frequencies in the range tested, which are quite different from the impedance 
of the test plate of very lightly damped structures. 
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Figure 25. Real part of tlie l(V' viscoelastic beam navegiiide absorber impedances 
at the center of the beam. 
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Figure 26. Iniaginary part of tlie 16" viscoelastic beam uaveguide absorber impe- 
dances at the center of beam. 
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Figure 27. Real part of the 16" constrained layer beam waveguide absorber impe- 
dances at the center of the beam. 
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IMPEDANCES VS. FREQUENCY 




Figure 28 . Imaginary part of the 16 " constrained layer beam naveguide absorber 
impedances at the center of beam. 
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Figure 29. Real part of the 20" constrained layer beam waveguide absorlier impe- 
dances at the center of the beam. 
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Figure 30. Imaginary part of tlie 20" constrained layer beam naveguide absorber 
impedances at the center of beam. 



C. DAMPING OF THE PLATE 

The purpose of this experiment is to investigate the dlect of light weiglit waveguide 
absorbers on damping increase in plate type structures over a wide frequency range. 
Damping of the test plate structure was measured with an attachment of each one of 
four different waveguide absorbers (16" or 20" viscoelastic or constrained layer beam) 
at three different locations. The frequency response and the modal damping values of the 
aluminum test plate with a waveguide absorber were compared with those of the same 
plate without any waveguide absorber over a wide frequency range (0-20U0Mz). 

The modal damping values of the test plate without a waveguide absorber are shown 
in Figure 31 on page 56 and Figure 32 on page 57. They show that more than 20 lightly 
damped modes are in the frequency range between lOOMz to 2000Mz with the lowest 
frequency 195IIz. 

The damping contribution of the same waveguide absorber to a same plate was dif- 
ferent depending on the location of the waveguide absorber attachment point. For the 
20" viscoelastic beam waveguide absorber, there was marked damping increase for 
modes lower than 1 1001 Iz when it was attached at location 1 (Figure 33 on page 58 and 
Figure 34 on page 59). However, it was more effective for modes higher than 1200Hz 
when it was attached at location 2 (Figure 35 on page 60 and Figure 36 on page 61). 
When it was attached at location 3 its contribution to plate damping increase was ef- 
fective for wider frequency range of 250Hz- 1500Hz (Figure 37 on page 62 and 
Figure 38 on page 63). This trend applies to other waveguide absorbers, though the 
magnitude of damping contribution and the frequency ranges were different for each 
case (Figures in Appendix B). 

Different modal damping effects with a different waveguide absorber at the same 
location were also observed. At location 1, the 20" viscoelastic beam waveguide absorber 
was most effective for the frequency range lower than llOOHz (Figure 33 on page 58 
and Figure 34 on page 59) and the 16" constrained layer beam waveguide absorber was 
most effective for the frequency range higher than 1200Hz ( Figure 39 on page 64 and 
Figure 40 on page 65). At location 2, the 20" and 16" viscoelastic beam waveguide ab- 
sorbers showed moderate contribution to plate damping and the 20" constrained layer 
beam waveguide absorber contributed the least amount to plate damping. I'or frequency 
ranges over l200FIz, the 16" constrained layer beam waveguide absorber pwoduced the 
largest damping increase. At location 3, the 20" and 16" viscoelastic beam waveguide 
absorbers and the 20" constrained layer beam waveguide absorber produced moderate 
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damping increase over a wide liequency range of 2dOHz- 1500Hz. However, the 16" 
constrained layer beam waveguide absorber showed the least damping increase. 

After the selection of optimum damping at each location, two wave guide absorbers 
(20" viscoelastic and 16" constrained layer beam) were attached at locations 1 and 2 and 
three waveguide absorbers (20" viscoelastic. 16" constrained layer and 16" viscoelastic 
beam) were attached at locations 1, 2 and 3 to see the effect of multiple waveguide ab- 
sorbers on the test plate damping (Figure 41 on page 66). Figure 42 on page 67 and 
Figure 43 on page 68 represent frequency respon.se and modal damping of the plate with 
two waveguide absorbers. Figure 44 on page 69 and Figure 45 on page 70 represent 
frequency response and modal damping of the plate with three waveguide absorbers. 
Damping of the plate increased with increasing number of the waveguide absorbers 
(Figure 42 on page 67 and Fugure 44 on page 69). The damping increase resulting from 
attachment of two or more absorbers resemble the sum of the seperatc result of each. 

Impedance of the test plate at location 1 is compared with the driving point impe- 
dance of the 20" viscoelastic beam in Figure 46 on page 71. This figure shows that at 
resonance frequencies the impedance of the test plate is much smaller than that of the 
waveguide absorber. Figure 47 on page 72 and Figure 48 on page 73 show the com- 
parsion between the impedances of the test plate at location 2 and location 3 with the 
drving point of the 16" constrained layer beam and 16" viscoelastic beam, respectively. 
These figures show that at low frequency range (<400Hz) the impedances of the test 
plate were higher than those of waveguide absorbers and at high frequency range 
(>400Hz) impedances of waveguide absorbers were higher at some resonance frequen- 
cies. 

I'rom equation (2.7), the loss factor contribution. ?/. of a waveguide absorber is 
proportional to the following 3 terms: 






|1 + 




f’ 



(4.1) 

(4.2) 



(4.3) 
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The first term is determined from waveguide characteristics only and tlie second term is 
determined from the impedance ratio between the test structure and a waveguide ab- 
sorber at the attachment point. However, tiie third term depends on the interaction be- 
tween the test structure and a waveguide absorber and the whole response of the test 
plate. IfZ^ is much smaller than Z., then will remain constant before and after the 
attachment of a waveguide absorber. Previous studies [Ref 2 and Ref 3J were based 
on this assumption and gave misleading indications that e.xperimental loss factors would 
be much higher than theoretical predictions. In the present study, as shown in 
Figure 46 on page 71. Figure 47 on page 72 and Figure 48 on page 73. the impedances 
of the test plate were similar in magnitude to impedances of w'aveguide absorbers at re- 
sonance frequencies and the loss factor contribution due to waveguide absorbers could 
not be predicted using equation (2.7). 
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Figure 31. The driving point frequency response of the test plate without a \>ave- 
guide absorber. 
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MODAL DAMPING VS. FREQUENCY 
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Figure 32 . 



Modal damping factors vs. frequency of the test plate without a wave- 
guide absorber. 
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FREQUBNCy (HZ) 




Figure 33. The driving point frequency response of the test plate vvitli a 20 " viscoe- 
lastic lieam waveguide ahsorher (dashed) at location 1 and without 
(solid). 
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MODAL DAMPING VS. FREQUENCY 
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Figure 34. Modal clamping factors vs. frequency of the test plate nitli a 2U" vis- 
coelastic beam waveguide absorber at location 1 and without. 
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FREQUENCY (UZ) 




Figure 35. The driving point frequency response of the test plate with a 20 " Aiscoe- 
lastic beam ■waveguide absorber (dashed) at location 2 and witliout 
(solid). 
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MODAL DAMPING VS. FREQUENCY 



Figme 36. 
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Modal damping factors vs. frequency of (he test plate with a 2U" viscoe- 
lastic beam waveguide absorber at location 2 and without. 
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FREQUENCY (HZ) 




(solid). 
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MODAL DAMPING VS. FREQUENCY 




Figure 38. Modal damping factors vs. frequency of the test plate ivitli a 20" viscoe- 
lastic beam waveguide absorber at location 3 and without. 
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FREQUENCY (HZ) 




without (solid). 
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Figure 40. Modal damping factors vs. frequency of the test plate with a 16" 
strained layer beam waveguide absorber at location 1 and without. 
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Figure 41. Tlie test plate uitli three waveguide absorbers. 
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waveguide absorber (dashed) and without (solid). 
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Figure 43. Modal damping factors vs. frequency of the test plate uitli a 20" viscoe- 
lastic (location 1) and a 16" constrained layer (location 2) beam wave- 
guide absorber and without. 
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lastic (at location 1), a 16" constrained layer (at location 2) and a 1(3' 
viscoelastic (at location 3) beam waveguide absorber (dashed) and 
without (solid). 
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Figure 45. Modal damping t'actors vs. frequency of tlie test plate nith a 20" viscoe- 
lastic (location 1), a 16" constrained layer (location 2) and a 16" viscoe- 
latic beam waveguide absorber and without. 
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V. COiNCLUSIONS 



The purpose of this study was to develop two kinds of high damping beam wave- 
guide absorbers for application to naval structures over a wide frequency range 
( 100Hz-2000Hz) and to evaluate theoretical prediction schemes using various beam 
theories. Comparison between theoretical predictions and experimental results show that 
Timoshenko beam theory and the six order beam theory can be used in the prediction 
of impedances of viscoelastic beam waveguide absorbers and constrained layer beam 
waveguide absorbers, respectively. 

Application of waveguide absorbers to a test plate structure showed that damping 
of the test structure can be increased significantly with one waveguide absorber over a 
frequency range. The magnitude and the frequency range of damping increase depend 
on the impedance of the waveguide absorber and on the location. The prediction of 
damping increase due to an attached waveguide absorber, equation (2.7). does not show 
the effect of waveguide absorber completely since this equation does not include the 
structure vibration reduction elfects due to the waveguide absorber. Since the contrib- 
ution of each waveguide absorber to the damping of the test plate is additive, the modal 
damping values of all the vibration modes of the test plate in a wide frequency range can 
be increased using multiple waveguide absorbers. 
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VL RECOMMENDATIONS 



From the results of present study follo’.ving are reconmiended for future studies on 
waveguide absorbers; 

1. Investigation of the eOect of a waveguide absorber on the reduction in vibration 
of the original structure and derivation of the relation between the impedances and the 
mode shape of the structure, 

2. Theoretical studies on impedance prediction of waveguides of dilferent shape such 
as circular viscoelastic plates and application to damping increase, 

3. Investigation of interaction between waveguide absorber for multiple waveguide 
absorbers applications, and 

4. Optimization of application of waveguide absorbers in design, number, and 
location. 
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APPENDIX A. 



\ 



1. THE IMPEDANCE OF THE INFINITE VISCOELASTIC BEAM USING 
BERNOULLI-EULER BEAM THEORY 



THE THREE REDUCE PROGRAMS WERE USED FOR EVALUATING THE IMPEDANCE 
OF THE INFINITE VISCOELASTIC BEAM USING BERNOULLI-EULER BEAM * 
THEORY. 



C 

C THE FIRST REDUCE PROGRAM DEFINED TRANSVERSE FORCE AND T\-}0 BOUNDARY 
C CONDITIONS. 

C 

YY: =Y1*EXP( I*K1*X)+Y3*E.'<P( I-'^K3*X); 

PP;=DF(YY,X); 

MM:=-EI*DF(PP,X); 

FF: =-DF(MM,X); 

Y0;=SUB(X=0,YY)-Y; 

P0:=SUB(X=0,PP)-P; 

BYE; 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Y1 AND Y3). 

C 

Y0;=Y1+Y3-Y; 

PO; =P1*Y1+P3*Y3-P; 

SOLVE(LST(YO,PO),Y1,Y3); 

BYE; 

C 

C THE THIRD REDUCE PROGRAM EVALUATED THE TRANSVERSE FORCE. 

C 

YY: =Y1*EXP( I*K1*X)+Y3*EXP( I*K3*X); 

FFF:=EI*DF(YY,X,3); 

FF:=SUB(X=0,FFF); 

BYE; 















THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE INFINITE 
VISCOELASTIC BEAM USING BERNOULLI-EULER BEAM THEORY. 



DEFINED 

E 

F 

FF 

G 

GRAV 

IMP 

POISS 

ROH 



VARIABLES; 

; COMPLEX YOUNG'S MODULUS OF THE VISCOELASTIC BEAM 
: FREQUENCY (HZ) 

: TRANSVERSE FORCE (LBF) 

: COMPLEX SHEAR MODULUS OF THE VISCOELASTIC BEAM 
: GRAVITY (33)6 IN/ SEC * 

: IMPEDANCE OF THE VISCOELASTIC BEAM 
: POISSON'S RATIO 

: DENSITY OF THE VISCOELASTIC BEAM 
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XB, XH,XL : DIMENSIONS OF THE BEAM (WIDTH, HEIGHT, LENGTH) 
XI : AREA MOMENT OF INERTIA 

XKl, XK2, XK3, XK4 : WAVE NUMBERS 

XNETA : ENERGY LOSS FACTOR OF THE VISCOELASTIC BEAM 
VI ; IMAGINARY (0,1) 









•;Wr 



c 

COMPLEX V1,G,E,KB1,KB2,XK1,XK2,XK3,XK4,P1,P3,A1,A2,A3,A4,FF,IMP 
REAL F,POISS,GRAV,PI,XB,XH,XG,ROH,XI 
REAL W , XLOG , XNETA , XA , XCA , XSA , KC , IMPR , IMPI 
C 

0PEN(UNIT=15,FILE=’IFI' ) 



CONSTANTS 



F=10. 0 
P0ISS=0.45 
GRAV=386. 0 
PI=3. 1415926 
V1=CMPLX(0. 0,1. 0) 

DIMENSIONS 



XB=1. 0 
XH=0. 35 



PROPERITIES 

R0H=0. 05505 /GRAV 
XI=(XB*XH**3)/12. 0 

FOUR DIFFERENT WAVE NUMBERS 

DO 100 1=1,100 
W=2*PI*F 

XL0G=ABS(L0G10(F/62)) 

XNETA=0. 65*EXP( -0. 52732*XL0G**1. 956) 

XNETA=0. 2 
XNETA=0. 5 

XG=0. 00002503*F**3-0. 1752*F**2+45 7. 5883*F+29280 
C XG=200000 

G=XG*CMPLX( 1. 0, XNETA) 

E=2*G*( 1+POISS) 

XA=ATAN( XNETA )/4. 0 
XCA=COS( -XA) 

XSA=SIN( -XA) 

KC=( SQRT( 1+XNETA**2) )**( -0. 25 ) 

KB 1=-CSQRT( ROH*XB*XH*W**2/ ( E*XI ) ) 

KB2= CSQRT( ROH*XB*XI W'*2/ ( E*XI ) ) 

XK1= CSQRT( KB 1 )*KC*CMPLX( XCA , XSA) 

XK2=-CSQRT( KB 1 )*KC*CMPLX( XCA , XS A ) 
XK3=-CSQRT(KB2)*KC*CMPLX(XCA,XSA) 

XK4= CSQRT( KB 2 ) *KC*CMPLX( XCA , XSA) 
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DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 



C 
C 

P1=V1*XK1 
P2=V1*XK2 
P3=V1-'>XK3 
P4=V1*XK4 
Al= P3/CP3-P1) 

A2=0 

A3=-P1/(P3-P1) 

A4=0 

FF= - V 1*E*X I* ( A 1*XK 1** 3+A 3*XK3 ** 3 ) 

C 

C IMPEDANCES 
C 

IMP=2*FF/(V1*W) 

IMPR=REAL(IMP) 

IMPI=AIMAG(IMP) 

WRITE(15,1000) F,IMPR,IMPI 
F=F+20. 0 
100 CONTINUE 

1000 F0RMAT(iX,D9. 3,2X,D12. 5,2X,D12. 5) 

CL0SE(UNIT=15) 

STOP 

END 

2 . THE IMPEDANCE OF THE FINITE VISCOELASTIC BEAM USING 
BERNOULLI-EULER BEAM THEORY 

YY: =YN^EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I->K3*X)+Y4*EXP( I-’-K4*X); 

PP: =DF(YY,X); 

MM: =-EI*DF(PP,X); 

FF: =-DF(MM,X); 

YO: =SUB(X=0,YY)-Y; 

P0:=SUB(X=0,PP)-P; 

ML: =SUB(X=L,MM); 

FL: =SU'3(X=L,FF)i 
BYE; 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Yl, Y2, Y3 
C AND Y4). 

C 

YO: =Y1+Y2+Y3+Y4-Y; 

PO: =P1*Y1+P2*Y2+P3*Y3+P4*Y4-P; 

ML: =MU'-Y1+M2*Y2+M3*Y3+M4*Y4; 

FL: =F 1*Y 1+F2*Y2+F3*Y3+F4*Y4; 

SOLVE ( LST( YO , PO , ML , FL) , Y 1 , Y2 , Y3 , Y4 ) ; 

BYE; 

C 

C THE THIRD REDUCE PROGRAM EVALUATED THE TRANSVERSE FORCE. 

C 

YY: =Y1*EXPC I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

FFF: =EI-'-DF(YY,X,3); 

FF: =SUB(X=0,FFF); 

BYE; 
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* 

* THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE FINITE 

* VISCOELASTIC BEAM USING BERNOULLI -EULER BEAM THEORY. 



7C 

Vf 



DEFINED VARIABLES; 



E 

F 



COMPLEX YOUNG'S MODULUS OF THE VISCOELASTIC BEAM 
FREQUENCY (HZ) 

TRANSVERSE FORCE (LfiFI 

COMPLEX SHEAR MODULUS OF THE VISCOELASTIC BEAM 
GRAVITY (32)5 IN/SEC 
IMPEDANCE OF THE VISCOELASTIC BEAM 
POISSON'S RATIO 

DENSITY OF THE VISCOELASTIC BEAM 
: DIMENSIONS OF THE BEAM (WIDTH, HEIGHT, LENGTH) 

AREA MOMENT OF INERTIA 
XK3, XK4 : WAVE NUMBERS 

ENERGY LOSS FACTOR OF THE VISCOELASTIC BEAM 
IMAGINARY (0,1) 

c 

COMPLEX V1,G,E,KB1,KB2,XK1,XK2,XK3,XK4,P1,P2,P3.P4 
COMPLEX M1,M2,M3,M4,F1,F2,F3,F4,DEN,A1,A2,A3,A4,FF,IMP 
REAL F , POISS , GRAV , P I , XB , XL , XH , XG , ROH , XI , W , XLOG , XNETA , XA , XCA , XSA 
REAL KC,KR1,KR2,KR3,KR4,KI1,KI2,KI3,KI4,IMPR,IMPI 
C 

0PEN(UNIT=15,FILE='FFI' ) 

CONSTANTS 



* 


FF : 


* 


G 


"jV 


GRAV 




IMP : 


•jV 


POISS : 


* 


ROH : 


•sV 


XE, XH,XL 


Vr 


XI ; 


* 


XKl, XK2, 


Vc 


XNETA 


* 


VI : 


* 





F=10. 0 
POISS=0. 45 
GRAV=386. 0 
PI=3. 1415926 
V1=CMPLX(0. 0,1. 0) 

DIMENSIONS 



XB=1. 0 
XL=8. 0 
XH=0. 35 



PROPERITIES 

ROH=0. 05 5 05 /GRAV 

FOUR DIFFERENT WAVE NUMBERS 

XI=(XB*.XH*---3)/12. 0 
DO 100 1=1,100 
W=2-’-PI*F 

XL0G=ABS(L0G10(F/62)) 

XNETA=0. 65*EXP( - 0. 52732*XLOG*-'l. 956) 
XNETA=0. 2 
XNETA=0. 5 
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XG=0. 00002503*F**3-0. 1752*F-v*2+45 7. 5883*F+29280 
C XG=200000 

G=XG*CMPLX( 1. 0 ,XNETA) 

E=2--'--G*(1+P0ISS) 

XA=ATAN(XNETA)/4. 0 
XCA=COS(-XA) 

XSA=SIN(-XA) 

KC=( SQRT( 1+XNETA’'*2) )**( -o. 25 ) 

KB 1=-CSQRT( R0H*XB*XH*W**2/ ( E*XI ) ) 

KB2= CSQRTC R0H*XB*XH*W**2/ ( E’>XI ) ) 

XK1= CSQRTC KB 1)*KC*CMPLX( XCA ,XSA) 
XK2=-CSQRT(KB1)*KC*CMPLX(XCA,XSA) 

XK3=-CSQRT( KB2)*KC*CMPLX( XCA , XSA) 

XK4= CSQRTC KB2)*KC*CMPLXC XCA , XSA) 

KR1=REALCXK1) 

KR2=REALCXK2) 

KR3=REALCXK3) 

KR4=REALCXK4) 

KI1=AIMAGCXK1) 

KI2=AIMAGCXK2) 

KI3=AIMAGCXK3) 

KI4=AIMAGCXK4) 

DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 



& 

& 

Ct 

& 

& 

& 

St 



P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

M 1=E XP C - XL*K 1 1 ) * C CO S C XL*KR 1 ) + V 1*S I NC XL*KR 1 ) ) *XK I** 2 
M2=EXP C -XL*KI 2 )*C COS C XL*KR2 ) +V1*S INC XL-'^KR2 ) )*XK2**2 
M3=EXPC -XL*KI 3)*C cose XL*KR3 )+Vl*SINC XL*KR3 ) )*XK3**2 
M4=EXPC-XL*KI4)*CCOSCXL*KR4)+Vl"SINCXL''KR4))*XK4*''-2 
F 1=EXP C -XL*K 1 1 ) * C COS C XL*KR 1 ) +V 1*S INC XL*KR 1 ) ) *XK1**3 
F2=EXP C -XL*KI 2 ) * C COS C XL’’- KR2 ) +V 1*S I NC XL*KR2 ) ) *XK2**3 
F3=EXPC -XL*K 1 3 ) * C COS C XL*KR3 ) +V I* S INC XL*KR3 ) ) *XK 3**3 
F4=EXPC-XL*KI4)*CC0SCXL*KR4)+V1*SINCXL*KR4))*XK4**3 
DEN=C F4*M3*P2-F4*M3*Pl-F4*M2*P3+F4*M2*Pl+F4*Ml*P3'F4*Mi*P2 
-F3*M4*P2+F3*M4*P1+F3*M2*P4-F3*M2*P1-F3*M1*P4+F3*M1*P2 



+F2*M4*P3-F2*M4*Pl-F2*M3*P4+F2*M3*Pl-‘-F2*Ml*P4-F2*Ml*P3 
-F1*M4*P3+F1*M4*P2+F1*M3*P4-F1*M3*P2-F1*M2*P4+F1*M2*P3) 
A1=C F4*M3*P2-F4*M2*P3-F3*M4*P2+F3*M2*P4't-F2*M4*P3-F2*M3*P4) 
/DEN 

A2=C-F4*M3*P1+F4*M1*P3+F3*M4*P1-F3*M1*P4-F1*M4*P3+F1*M3*P4) 

/DEN 



A3=C F4*M2*P 1 -F4*M 1*P2 -F2*M4*P 1+F2*M 1*P4+F 1*M4*P2 -F1*M2*P4 ) 
/DEN 

A4=C -F3*M2*P 1+F3*M1*P2+F2*M3*P 1 -F2*M 1*P3 -F l*M3*P2-f-F 1*M2*P3 ) 
/DEN 

FF=-V1*E*XI*CA1*XK1**3+A2*XK2**3+A3*XK3**3+A4*XK4**3) 



IMPEDANCES 



IMP=2*FF/CV1*W) 

IMPR=REALCIMP) 

IMPI=AIMAGCIMP) 



SO 



WRITE( 15,1000) F,IMPR,IMPI 
F=F+20. 0 
100 CONTINUE 

1000 FORMAK IX, D9. 3 , 2X, D12. 5 , 2X,D12. 5) 

CL0SE(UNIT=15) 

STOP 

END 

3. THE IMPEDANCE OF THE FINITE VISCOELASTIC BEAM USING 



TIMOSHENKO BIiAM THEORY 



"jV V? V> Vr Vr V * 'k ir Vr V? Vr k 



* V * V Vc * V k *3 V k “V 



■sWriWcVc” 



•;W?*jV*,V*/Wc’ 



VftWr'jWc 






* THE THREE REDUCE PROGRAMS WERE USED FOR EVALUATING THE IMPEDANCE 

* OF THE FINITE VISCOELASTIC BEAM USING TIMOSHENKO BEAM THEORY. * 

•jV *V 






V r nV nV iV “.V "“.V V? r.V “V ^ I- r 



r :V • ' k / f "j » •} V *j ■ 






C 

C THE FIRST REDUCE PROGRAM DEFINED SHEAR FORCE AND FOUR BOUNDARY 
C CONDITIONS. 

C 

YY: =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

PPl: =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

PP2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X); 

PP;=PP1+PP2; 

MM: =DF(PP,X); 

SS:=-(DF(YY,X)-PP); 

YO: =SUB(X=0,YY)-Y; 

P0:=SUB(X=0,PP)-P; 

ML: =SUB(X=L,MM); 

SL:=SUB(X=L,SS); 

BYE; 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Yl, Y2, 

C Y3 AND Y4). 

C 

YO: =Y1+Y2+Y3+Y4-Y; 

PO : =R 1*Y 1 +R2*Y2+R3*Y3+R4* Y4 - P; 

ML: =M1*Y 1+M2*Y2+M3*Y3+M4*Y4; 

SL: =S 1*Y1+S2*Y2+S3*Y3+S4*Y4; 

SOLVE(LST(YO,PO,ML,SL) ,Y1,Y2,Y3,Y4); 

BYE; 

C 

C THE THIRD REDUCE PROGRAM EVALUATED THE SHEAR FORCE. 

C 

YY: =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

PPl: =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

PP2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4>’-X); 

PP:=PP1+PP2; 

S S S : = - KA*GG* AA* ( DF ( Y Y , X ) - PP ) ; 

SS:=SUB(X=0,SSS); 

BYE; 



* * 
* THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE FINITE * 
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VISCOELASTIC BEAM USING TIMOSHENKO BEAM THEORY. 



* 



DEFINED VARIABLES; 

COMPLEX YOUNG'S MODULUS OF THE VISCOELASTIC BEAM 
FREQUENCY (HZ) 

COMPLEX SHEAR MODULUS OF THE VISCOELASTIC BEAM 
GRAVITY (32^ IN/ SEC 
IMPEDANCE OF THE VISCOELASTIC BEAM 
SHAPE OF CROSS SECTION 
POISSON'S RATIO 

DENSITY OF THE VISCOELASTIC BEAM 
SHEAR FORCE (LBF) 

XB, XH,XL : DIMENSIONS OF THE BEAM (WIDTH, HEIGHT, LENGTH) 
XI : AREA MOMENT OF INERTIA 

XKl, XK2, XK3, XK4 : WAVE NUMBERS 

XNETA : ENERGY LOSS FACTOR OF THE VISCOELASTIC BEAM 
VI : IMAGINARY (0,1) 



E 

F 

G 

GRAV 

IMP 

KAR 

POISS 

ROH 

SS 



Vf **3WrV? Vr ‘V* 









c 

COMPLEX V 1 , G , E , XKNl , XKN2 , XKD , XK 1 , XK2 , XK3 , XK4 , R1 , R2 , R3 , R4 
COMPLEX Ml , M2 , M3 ,M4 , S 1 , S2 , S3 , S4 ,DEN . A1 , A2 , A3 , A4 , SS , IMP 
REAL F , KAR , PO I S S , GRAV , P I , XB , XL , XH , XG , ROH , X I , W , XLOG , XNETA 
REAL KRl , KR2 , KR3 , KR4 , KI 1 , KI 2 , KI 3 , KI4 , IMPR , IMPI 
C 

OPEN(UNIT=15,FILE='TIM' ) 

CONSTANTS 

F=10. 0 
KAR=0. 83 
POISS=0. 45 
GRAV=386. 0 
PI=3. 1415926 
V1=CMPLX(0. 0,1. 0) 

DIMENSIONS 



XB=1. 0 
XL=8. 0 
XH=0. 35 



PROPERITIES 

R0H=0. 05505/GRAV 
XI=(XB*XH*-->3)/12. 0 

FOUR DIFFERENT WAVE NUMBERS 

DO 100 1=1,100 
W=2-'-PI*F 

XLOG=ABS(LOG10(F/62)) 

XNETA=0. 65"EXP( -0. 52732*XLOG**l. 956) 
XNETA=0. 2 
XNETA=0. 5 
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XG=0. 00002503*F**3-0. 1 75 2* F''* 2+45 7. 5883*F+29280 
C XG=200000 

G=XG*CMPLX( 1. 0 ,XNETA) 

E=2+G*(1+P0ISS) 

XKNl=W**2-'-( KAR-'-G+E ) /ROH 
XKN2= SQRT( W**4*( ( KAR*G+E) /R0H)**2 

-4*KA R*G*E* ( W**4 - KAR*G*XB->XH*W** 2 / ( ROH*X I ) ) 

/(ROH**2)) 

XKD=2*KAR*G*E/( ROH**2 ) 

XK1= CSQRT((XKN1+XKN2)/XKD) 

XK 2= - C SQRT ( ( XKN 1+XKN 2 ) / XKD ) 

XK3= CSQRT((XKN1-XKN2)/XKD) 

XK4=-CSQRT((XKN1-XKN2)/XKD) 

KR1=REAL(XK1) 

KR2=REAL(XK2) 

KR3=REAL(XK3) 

KR4=REAL(XK4) 

KI1=AIMAG(XK1) 

KI2=AIMAG(XK2) 

KI3=AIMAG(XK3) 

KI4=AIMAG(XK4) 

DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 

R 1=( ROH*W** 2 -KAR*G*XK 1**2 ) / ( V 1*KAR*G*XK 1 ) 

R2=( ROH*W**2 -KAR*G*XK2**2 ) / ( V 1*KAR*G*XK2 ) 

R3=( ROH*W**2-KAR*G*XK3**2 ) / ( V 1*KAR*G*XK3 ) 

R4=( ROH*W**2 -KAR*G*XK4** 2 ) / ( V1*KAR*G*XK4 1 
M1=EXP( -XL*KI 1 )*( COS( XL*KiU )+Vl*SIN( XL*KR1 ) )*V1*XK1*R1 
M2=EXP( -XL*KI2 )*( COS( XL*KR2 )+Vl*SIN( XL*KR2 ) )*V1 *XK2*R2 
M3=EXP( -XL*KI3 )*( COS(XL*KR3 ) +V1*SIN( XL*KR3) )*V1*XK3*R3 
M4=EXP ( -XL*K 1 4 ) * ( C 0 S ( XL*KR4 ) +V 1* S I N ( XL*KR4 ) ) * V 1*XK4*R4 
S1=EXP( -XL*KI l)*(COS( XL*KR1 )+Vl*SIN(XL*KRl) )*( R1 -V1*XK1 ) 
S2=EXP( -XL*KI 2)*( COS( XL*KR2 )+Vl*S IN( XL*KR2 ) )*( R2 -V1*XK2 ) 
S3=EXP( -XL*KI 3 ) *( COS( XL*KR3 ) +V1*S IN( XL*KR3 ) ) *( R3 - V 1*XK3 ) 
S4=EXP(-XL*KI4)*(C0S(XL*KR4)+V1*SIN(XL*KR4))*(R4-V1*XK4) 

DEN=( S4*M3*R2 - S4*M3*R 1 - S4*M2*R3+S4*M2*R 1+S4*M 1*R3 - S4*M 1*R2 
- S3*M4*R2+S3*M4*R 1+S3*M2*R4- S3*M2*R 1 - S3*M1*R4+S3*M 1*R2 
+S2*M4*R3-S2*M4*R1-S2*M3*R4+S2*M3*R1+S2*M1*R4-S2*M1*R3 
- S 1*M4*R3+S 1*M4*R2+S 1*M3*R4 -S 1*M3*R2 - S 1*M2*R4+S 1*M2*R3 ) 

A 1=( S4*M3*R 2 - S4*M2*R3 - S 3*M4*R2+S 3*M2*R4+S 2*M4*R 3 - S 2*M3*R4 ) 

& /DEN 

A2=( - S4*M3*R 1+S4*M 1*R3+S 3*M4*R 1 - S3*M1*R4 - S 1*M4*R3+S 1*M3*R4 ) 

& /DEN 

A3=( S4*M2*R 1 - S4*M 1*R2 - S 2*M4*R 1+S2*M 1*R4+S 1*M4*R2 - S 1*M2*R4 ) 

& /DEN 

A4=( -S3*M2*R1+S3*M1*R2+S2*M3*R 1 -S 2*M1*R3 -S 1*M3*R2+S 1*M2*R3 ) 

8c /DEN 

SS=KAR*G*XB*XH*( Al*( R1 -V1*XK1 )+A2*( R2-V1*XK2 )+A3*( R3-V1*XK3 ) 
6c +A4*(R4-V1*XK4) ) 

C 

G IMPEDANCES 
C 

IMP=2*SS/(V1*W) 

IMPR=REAL(IMP) 

IMPI=AIMAG(IMP) 
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WRITE(15,1000) F,IMPR,IiMPI 
F=F+20. 0 
100 CONTINUE 

1000 F0RMAT(1X,D9. 3 , 2X ,D12. 5 , 2X,D 12. 5) 

CLOSE(UNIT=15) 

STOP 

END 

4. THE IMPEDANCE OF THE FINITE CONSTRAINED LAVER BEAM 
USING THE SIXTH ORDER BEAM THEORY 



THE THREE REDUCE PROGRAMS WERE USED FOR EVALUATING THE IMPEDANCE * 
OF THE INFINITE VISCOELASTIC BEAM USING BERNOULLI -EULER BEAM 
THEORY. 



C 

C THE FIRST REDUCE PROGRAM DEFINED TRANSVERSE FORCE AND 'n-ZO BOUNDARY 
C CONDITIONS. 

C 

YYl: =Y1*EXP( I*K1*X)+Y2*EXP( I*K2->X); 

YY2: =Y3*EXP( I*K3-’-X)+Y4*EXP( I*K4*X); 

YY3; =Y5*EXP( I*K5^^X)+Y6*EXP( I*K6*X); 

UUl: =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

UU2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X); 

UU3: =R5*Y5--’^EXP( I*?;5*X)+R6*Y6*EXP( I*K6*X); 

YYY:=YY1+YY2+YY3; 

UUU: =UU1+UU2+UU3; 

YO: =SUB(X=0,YYY)-Y; 

P0:=SUB(X=0,DF(YYY,X))-Pi 
UO: =SUB(X=0,UUU)-U; 

S L; =SUB ( X=L , ( PP*DF ( YYY , X . 3 ) +QQ*DF( YYY , X ) +RR*UUU ) ) ; 

ML: =SUB ( X=L , ( PP*DF( YYY , X , 2 ) +SS*DF( UUU , X) ) ) ; 

NL: =SUB( X=L , ( TT*DF( UUU ,X) ) ) ; 

BYE; 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Y4, Y5 , 

C AND Y6). 

C 

YO: =Y1+Y2+Y3+Y4+Y5+Y6-Y; 

P0:=P1*Y1+P2*Y2+P3*Y3+P4*Y4+P5*Y5+P6*Y6-P; 

UO: =R1->Y1+R2*Y2+R3*Y3+R4*Y4+R5*Y5+R6’'-Y6-U; 

SL:=S1*Y1+S2*Y2+S3*Y3+S4*Y4+S5*Y5+S6='-Y6; 

^lL:=Ml*Yl+M2■'•■Y2+M3•■•Y3+^14•'•-Y4+M5’’•Y3-^•M6*Y6; 

NL: =N U'- Y l+N2-*Y2+N3-'^Y3+N4*Y4+N5*Y5FN6-'-Y6; 

SOLVE (LST(SL, ML, NL) ,Y4,Y5 ,Y6); 

BYE; 

C 

C THE THIRD REDUCE PROGRAM SOLVED VARIABLES (Yl, Y2 AND Y3) USING 

G Y4, Y5 AND Y6. THE THIRD REDUCE PROGRAM EVAJ.UATED THE SHEAR FORCE. 

C 

Y4: =BU>Y1+C1*Y2+D1*Y3; 

Y5: =B2*Y1+C2*Y2+D2*Y3; 

Y6: =B3*Y1+C3*Y2+D3*Y3; 

YO: =Y1+Y2+Y3+Y4+Y5+Y6-Y; 
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PO: =P1*Y1+P2*Y2+P3*Y3+F4*Y4+P5*Y5+P6*Y6-P; 
UO: =R1*Y1+R2*Y2+R3*Y3+R4“ Y4+R5-'=-Y5+R6'-Y6 -U; 
SOLVE(LST(YO,PO,UO) ,Y1,Y2,Y3); 

YYYl: =Y1--'EXF( I*K1->X)+Y2*EXP( I*K2*X); 

YYY2: =Y3*EXP( I'>K3*X)4-Y4*EXP( I*K4*X) ; 

YYY3: =Y5*EXP( I*K5*X)4-Y6-'«EXP( I*K6*X) ; 

YYYY: =YYY1+YYY2+YYY3; 

UUUl: =R1*Y1*EXP( I’>K1*X)H-R2*Y2*EXP( I*K2*X) ; 
UUU2: =R3'-Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X) ; 
UUU3: =R5*Y5*EXP( I*K5*X)+R6*Y6*EXP( I*K6*X) ; 
UUUU: =UUU1+UUU2+UUU3; 

SSS: =PP*DF( YYYY , X , 3 ) +QQ*DF( YYYY , X)+RR*UUUU; 
SS:=SUB(X=0,SSS); 

BYE; 



•V vV 

* THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE CONSTRAINED * 

* LAYER BEAM USING THE 6TH ORDER BEAM THEORY. 



•)C 

* 

* 



Vc 

* 

Vf 

Vc 

* 

* 

* 

* 

Vf 

•jV 



DEFINED VARIABLES; 

El, E3 : COMPLEX YOUNG'S MODULUS OF THE CONSTRAINED LAYER 

F : FREQUENCY (HZ) 

G : COMPLEX SHEAR MODULUS OF THE VISCOELASTIC 

GRAV : GRAVITY (386 IN/SEC**2) 

HI, H3 : HEIGHT OF THE CONSTRAINED LAYER 

H2 : HEIGHT OF THE VISCOELASTIC 

II, 12 ; AREA MOMENT OF INERTIA 

IMP ; IMPEDANCE OF THE CONSTRAINED LAYER BEAM 

KAR : SHAPE OF CROSS SECTION 

POISS : POISSON'S RATIO 

ROHl : DENSITY OF THE CONSTRAINED LAYER 

ROH2 : DENSITY OF THE VISCOELASTIC 

SS : SHEAR FORCE (LBF) 

XB, XL : DIMENSIONS OF THE BEAM (WIDTH, LENGTH) 

XKl, XK2, XK3, XK4, XK5 , XK6 : WAVE NUMBERS 

XNETA : ENERGY LOSS FACTOR OF THE VISCOELASTIC BEAM 

VI : IMAGINARY (0,1) 



if 



c 

COMPLEX*32 V1,G,BB,P1,P2,P3,P4,P5,P6,R1,R2,R3,R4,R5,R6,QQ,RR 
COMPLEX-'-32 S 1 , S2 , S3 , S4 , S5 , S6 , Ml , M2 , M3 , M4 , M5 , M6 , N1 , N2 , N3 , N4 , N5 , N6 
COMPLEX*32 B1,C1,D1,B2,C2,D2,B3,C3,D3,DEN1,DEN2,A1,A2,A3 
COMPLEX B , D , XKl , XK2 , XK3 , XK4 , XK5 , XK6 . XXKl , XXK2 , XXK3 , IMP 
REAL*16 F,PI,GRAV,XB,XL,XG,H1,H2,H3,E1,E3,I1,I3,DD,R0H1,R0H2 
REAL*16 CC,XM,T1(30) ,T2(30) ,T3(30) ,Q1(30) ,Q2(30) ,Q3(30) 

REAL*16 XL0G,W,Z,C,BR,BI,CR,DR,DI,FT1,FQ1,D1T1,DTQ1,DQQ1,UQT1 
REAL* 16 FT2 , FQ2 , DTT2 , DTQ2 , DQQ2 , DQT2 , FT3 , FQ3 , DTT3 , DTQ3 , DQQ3 , DQT3 
REAL* 16 XKR 1 , XKR2 , XKR3 , XKI 1 , XKI 2 , XKI 3 , KRl , KR2 , KR3 , KR4 , KR5 , KR6 
REAL*16 KI1,KI2,KI3,KI4,KI5,KI6,PP,S3,TT,IMPR,IMPI 
REAL XNETA 
C 

0PEN(UNIT=15,FILE='SAN' ) 

C 
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CONSTANTS 



F=2010 

PI=3. 1415926 
GRAV=386. 0 
V1=CMPLX(0. 0,1. 0) 

DIMENSIONS 



XB=1. 0 
XL=8. 0 
H1=0. 0625 
H2=0. 065 
H3=0. 0625 

PROPERITIES 

El=10000000. 0 
£3=10000000. 0 
I1=XB-->H1**3/12. 

I3=XB->H3-”-*3/12. 

DD=E1*I1+E3*I3 
R0H1=0. 09832/GRAV 
R0H2=0. 03663/GRAV 
CC=Hl/2. +H2+H3/2. 

XM=XB*( ( H1+H3 ) •'■-R0H1+H2*R0H2 ) 

INITIAL GUESS WAVE NUMBERS 

Tl( l)=-4. 31911862 
T2( l)=-0. 05139603 
T3( 1)= 3. 77697017 
Ql(l) = -0. 33835599 
Q2( l)=-0. 05104052 
Q3(l)=-0. 27669665 

FIND SIX DIFFERENT WAVE NUMBERS USING NEWTON'S METHOD 
SIX DIFFERENT WAVE NUMBERS ARE XKl ,XK2 ,XK3 ,XK4 ,XK5 ,XK6 

DO 100 1=1,100 
W=2^’^PI-'-F 

XLOG=ABS(LOG10(F/800)) 

XNETA=1. 1*EXP( -0. 3906*XL0G->*1. 53) 

XNETA=0. 5 
XNETA=1. 0 

XG=-0. 00001258*F**2+0. 2472*F+74. 988 
C XG=1000 

G=XG*CMPLX( 1. 0 ,XNETA) 

B B=G*( E 1-’^H1+E3*H3 ) / ( H2*E l'''H l^>E3*iI3 ) 

Z=CC**2*E l^’^H 1*E 3^'iI3 / ( DD* ( E l^’-H 1+E 3 ■•-}I3 ) ) 

B=BB*(1. 0+Z) 

C=-XM^'^W**2/DD 

D=-BB*XM*W**2/DD 

BR=REAL(B) 

BI=AIMAG(B) 

CR=C 
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DR=REAL(D) 

DI=AIMAG('D) 

DO 200 J=l,15 

FT 1=T 1 ( J ) ** 3+B R*T 1 ( J ) ** 2 - 3*Q 1 ( J ) * -’'2 - 2‘>B I *Q 1 ( J ) +C R ) 1 ( J ) 

& +(-BR*QUJ)**2+DR) 

FQ 1=Q 1 ( J ) ** 3+B I*Q 1 ( J ) *-->2+( - 3*T 1 ( J ) ** 2 - 2 ->BR*T 1 ( J ) - CR ) '->0 1 ( J ) 
5c +(-BI*Tl( J)**2-UI) 

DTT 1=3*T 1 ( J )■•■■■•■■ 2+2*B R* T 1 ( J ) 

5c +( -3*Q1( J)*---2-2’'-BI*Ql( J)+CR) 

DTQ 1=( - 6*T 1 ( J ) - 2*BR ) *Q 1 ( J ) - 2*B I*T 1 ( J ) 

DQQ1=3*Q1( J)**2+2*BI*Q1(' J) 

5c +(-3*Tl(J)**2-2*BR*Tl( J)-CR) 

DQT 1=( - 6 *Q 1 ( J ) - 2* B I V^Tl ( J ■) - 2* B R*0 1 ( J ■) 

T 1( J+ 1 ) =T 1 ( J ) - ( DQQ 1*FT 1 -DTQ I* FQ 1 ) / ( DTT 1*DQQ 1 -DTO 1 -'-DOT 1 ) 

Q 1 ( J+ 1 ) =Q 1 ( J ) - ( -DQT l*FT i+DTT 1*FQ 1 ) / ( DTT 1*DQQ 1 - DTO 1*D0T1 ) 
FT2=T2 ( J) 3+BR*T2 ( J ) *-'-2+( - 3 '■•Q2 ( J ) **2 - 2*B I >'-Q2 ( J ) +GR ) '>T2 ( J ) 
5c +(-BR*Q2(J)**2+DR) 

FQ2=Q2( J)**3+BI*Q2( J)*''-2+f-3*T2(J)**2-2*BR-'-T2(J)-CR)*Q2(J) 
5c +(-BI*T2(J)**2-DI) 

Dtt2=3*t2( J)*’V2+2*BR*T2( J) 

Sc +(-3*02(J)**2-2*BI*Q2(J)+CR) 

DTQ2=( - 6*T2 ( J ) - 2*B R ) Q2 ( J ) - 2*B I*T2 ( J ) 

DQQ2=3*Q2( J)**2+2*BI*Q2( J) 

5c +( -3*T2( J)**2- 2*BR'’-T2( J) -CR) 

DQT2=( - 6*Q2 ( J ) - 2* B I ) *T2 ( J ) - 2*BR*Q2 ( J ) 

T2( J-l-1 )=T2( J) - ( DQQ2*FT2 -DT02*FQ2 ) / ( DTT2*DQQ2-DTQ2*D0T2 ) 
Q2( J+1 ) =Q2 ( J) - ( -DQT2---FT2+DTT2*FQ2 ) / ( DTT2*DQQ2 -DT02*DQT2 ) 
FT3=T3 ( J ) ** 3+BR*T3 ( J) 2+( - 3 ->03 ( J ) **2 - 2*B I ■■'Q3 ( J ) +CR) *T3 ( J ) 
5c +(-BR*Q3(J)**2+DR) 

FQ3=Q3( J)**3+BI*Q3( J)**2+( -3-'-T3( J)**2-2’''BR*T3( J) -CR)-'-Q3( J) 
5c +( -BI*T3(J)**2-DI) 

DTT3=3*T3( J) **2+2*BR*T3( J) 

5c +( -3*Q3( J)**2-2*BI*Q3( J)+CR) 

DTQ3=( -6*T3( J) -2*BR)*Q3( J) -2*BI*T3( J) 

DQQ3=3*Q3 ( J)**2+2*B I-’-Q3 (, J) 

5c +( -3*T3( J)**2-2*BR'''T3( J) -CR) 

DQT3=( - 6*Q3 ( J ) - 2*B I )*T3 ( J ) - 2*BR*Q3 ( J ) 

T3( J+1 )=T3( J) - ( DQQ3*FT3-DTQ3*FQ3 ) / ( DTT3*DQQ3 -DTQ3*DQT3 ) 
Q3( J+1 )=Q3( J) - ( -DQT3*FT3+DTT3*FQ3 ) / ( DTT3*DQQ3 -DTQ3*DQT3 ) 
200 CONTINUE 

XKR1= T1(J) 

XKR2= T2(J) 

XKR3= T3(J) 

XKI1= Q1(J) 

XKI2= Q2(J) 

XKI3= Q3(J) 

XXK 1=CMPLX( XKR 1 , XKI 1 ) 

XXK2=CI1PLX( XKR2 , XKI2 ) 

XXK3=CMPLX(XKR3,XKI3) 

XK1= CSQRT(XXKl) 

XK2=-CSQRT(XXK1) 

XK3= CSQRT(XXK2) 

XK4=-CSQRT(XXK2) 

XK5= CSQRT(XXK3) 

XK6=-CSQRT(XXK3) 

KR1= REAL(XKl) 
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KE2= REAL(XK2) 

KR3= REAL(XK3) 

KR4= REAL(XK4) 

KR5= REAL(XK5) 

KR6= REAL(XK6) 

KI1=AIMAG(XK1) 

KI2=AIMAG(XK2) 

KI3=AIMAG(XK3) 

KI4=AIMAG(XK4) 

KI5=AIMAG(XK5) 

KI6=AIf!AG(XK6) 

DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 

P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

P5=V1*XK5 

P6=V1*XK6 

Rl=Vl*BB*Z*DD*XKl/( ( XK1**2+BB)*E3*H3*CC) 

R2=V 1*B B*2*DD*XK2 / ( ( XK2 -’'*2+B B ) *E 3*H3 -’'CC ) 
R3=Vl*BB*Z-'-DD*XK3/((XK3**2+BB)*E3*H3*CC) 

R4=V 1*BB*Z*DD*XK4 / ( ( XK4’>*2+BB )*E3*H3*CC ) 
R5=Vl-''BB*Z*DD*XK5/( ( XK5**2+BB )*E3*H3*CC ) 
R6=Vl*BB*Z-’-DD*XK6/( ( XK6->*2+BB)*E3*H3*CC-) 

PP=DD 

QQ=-DD*BB-'-Z 

RR=DD*BB*Z*( E 1*H1+E3*H3 ) / ( CC-'-E 1*H1 ) 

SS=-CC*E3*H3 

TT=E3*H3 

S 1=( -V1*PP*XK1**3+V1*QQ*XK1+RR*R1 )*EXP( -XL*KI 1 ) 

& *(C0S(XL*KR1)+V1*SIN(XL*KR1) ) 

S2=( -Vl*PP'’-XK2**3+Vl*QQ-’-XK2H-RR*R2)*EXP( -XL*KI2) 

&. *( COS( XL*KR2 )+Vl*SIN( XL^'-KR2) ) 

S3=( -V 1*PP*XK3**3+V l’'-QQ^''XK34-RR*R3 ) *EXP( -XL*KI 3 ) 

5c *( COS( XL*KR3 )+Vl*S IN( XL*KR3 ) ) 

S4=( -Vl-’--PP-'-XK4**3+Vl*QQ*XK4+RR*R4)-'^EXP(-XL-'^KI4) 

Sc *( C0S(XL*KR4)+V1->SIN(XL->KR4) ) 

S5=( -Vl*PP*XK5*--^3+Vl*QQ->XK5+RR^’-R5 ) *EXP( -XL^-'KI5 ) 

Sc *(C0S(XL^'^KR5)+V1*SIN(XL*KR5) ) 

S6=( -Vl^’^PP*XK6-->*3+Vl*0Q*XK6+RR'-R6)’'^EXP( -XL-’^KI6) 

& *( COS( XL*KR6)+V1*SIN( XL*KR6) ) 

Ml=( V1*SS*R1*XK1 -PP*XK1*^'' 2)-”-EXP( -XL*KI 1 ) 

Sc * ( COS ( XL*KR 1 ) +V 1* S I N ( XL*KR 1 ) ) 

M2=( V1’^SS*R2*XK2-PP*XK2**2)*EXP( -XL*KI2 ) 

Sc *(COS(XL^'-KR2)+Vl*SIN(XL*KR2) ) 

M3=( V 1*SS*R3*XK3 - PP*XK3*'''2 J*EXP( -XL*K 1 3 ) 

Sc *( C0S( XL*KR3 ) +V 1*S IN( XL*KR3 ) ) 

M4=( V1*SS*R4*XK4-PP*XK4**2)*EXP( -XL*KI4) 

Sc *(COS( XL’’-KR4)+V1*SIN( XL*KR4) ) 
M5=(V1’'-SS*R5*XK5-PP*XK5**2)-"'EXP(-XL’>KI5) 

Sc * ( C OS ( XD>KR5 ) +V 1* S J. N ( XL*KR5 ) ) 

M6=( V1*SS*R6^'--XK6 -PP’'-'XK6**2 )*EXP( -XL*KI 6 ) 

Sc *( COS ( XL*KR6 ) +V 1*S IN( XL*KR6 ) ) 

N 1=V 1*TT*R 1*XK l^'- EXP ( - XL*K 1 1 ) 
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p^p'p'p^p^p’p^p^p'p’p^ p’p^p’p’p’p^p'p’ 



& *( cos ( XL*KR 1 ) +V 1*S IN( XL^>KR1 ) ) 

N2=V1*TT*R2*XK2*EXP( -XL*KI 2 ) 

& *( C0S(XL*KR2)+V1*SIN(XL*KR2) ) 

N3=Vl*Tr-'R3*XK3*EXF(-XL^’^KI3) 

& 'nC0S(XL*KR3)+Vl*SIN(XL*KR3) ) 

N4=V 1 -■'TT-'R4*XK4*EXP ( - XL*KI 4 ) 

& *(C0S(XL*KR4)+V1*SIN(XL*KR4) ) 

N5=V1’^TT*R5*XK5-->EXP( -XL*KI5 ) 

& *(C0S(XL*KR5)+V1*SIN(XL*KR5) ) 

N6=V1*TT*R6*XK6*EXP( -XL*KI 6 ) 

& *( COS ( XL*KR6 ) +V 1*S IN( XL-'-KR6 ) ) 

Bl=( -N6*M5'>S1+N6*M1*S5+N5*M6^'^S1-N5*M1*S6-N1*M6*S5+N1*M5*S6) 
& /(N6*M5-''S4-N6-'-M4'VS5-N5*M6*S4+N5*M4*S6+N4*M6'''S5-M4*M5*S6) 

Cl=( -N6*M5*S2+N6'-M2*S5fN5*M6*S2-H5*M2*S6-N2’-M6*S5+N2-'M5-'S6) 
& / ( N6*M5*S4 -N6*M4*S5 -N5*M6*S4+N5*iM4*S6+N4->M6*S5 -N4*I15-''S6 ) 

Dl=( -N6*M5*S3+H6*M3*S5+N5*M6*S3-N5*M3*S6-N3*M6-'-S5+N3-’-M5*S6) 
& /(N6*M5*S4-N6*M4*S5-N5'’-M6*S4+N5*M4-'-S6+N4-’-M6>'-S5-N4-’-M5'’-S6) 

B2=( N6*M4*S 1 -N6*M1*S4-N4*M6*S 1+N4*M1*S6+N1*M6*S4-N1*M4>>S6 ) 

& / ( N6*M5*S4-N6*M4*S5 -N5-'-M6’''S4+N5*M4*S6+N4*M6*S5 -N4*H5*S6 ) 

C2=( N6*M4* S 2 -N6*M2*S4 -N4*M6*S 2+M4*M2=>S6+N2-’^M6 '>S4 -N2*M4*S 6 ) 

& / ( N6*M5*S4-N6*M4*S5 -N5*M6*S4+N5*M4*S6+N4*tl6*S5 -N4->M5*S6 ) 

D2=( N6-'-M4*S3 -N6*M3*S4 -N4*M6*S3+N4*M3*S6+N3*M6*S4 -N3*M4*S6 ) 

& /(N6*M5*S4-N6*M4*S5-N5*tl6*S4+N5*M4*S6+N4*M6*S5-N4*M5*S6) 

B3=( -N5’>M4*Sl+N5*Ml*S4+N4*M5*Sl-N4*Ml--'-S5-Nl-'-M5*S4+Nl-->fl4*S5) 
& /(N6*M5*S4-N6*M4’>S5-N5*M6*S4+N5*M4*S6+N4*M6*S5-N4*M5’'S6) 

C3=( -N5*M4*S2+N5*M2*S4+N4’’-M5--''S2 -N4*M2*S5 -N2''-M5’-S4+N2*M4*S5 ) 
& / ( N6*M5*S4 -Ne'^W^SS -N5*tl6*S4+N5*M4*S6+N4*M6*S5 -N4*M5*S6 ) 

D3=( -N5*M4*S3+N5*t13*S4+N4*M5*S3-N4*M3*S5-N3*M5*S4+N3*M4*S5) 
& /(N6*M5*S4-N6*M4*S5-N5*M6’>S4+N5*M4*S6+N4*M6*S5-N4*M5*S6) 

DEN1=( R6*P5 -R5*P6)*( D3*C2*B 1+D3*C2-D3*B2*C1-D3*B2-C3*D2*B1 

-C3*D2+C3*B2'’-D 1+C3*B2+B3*D2*C 1+B3*D2 
-B3^>C2'’'Dl-B3-’-C2) 

+( R4*P6 -R6*P4 )*( D3*C2*B 1 -D3''-B2*C 1 -D3*C 1+D3*B1 -C3’’'D2*B 1 
+C3*B 2 •• D 1+C 3*D 1 - C3*B l+B 3*D2*C 1 - B 3->C2 -'D 1 
-B3"D1+B3“C1) 

+( R5*P4-R4*P5 ) *( D3-’'C2^'-B 1 -D3*B2*C1 -C3*D2*B 1+C3*B2*D1 
+B3*D2*C 1 -B3*C2*D] +D2*C 1 -D2*B 1 -C2*D 1 
+C2-’'B1+B2*D1-B2*C1) 

DEN2=( R3*P6 -R6*P3 ) *( C3*B 2+C3*B 1+C3 -B 3*C2 -B 3'->C 1 -B3 ) 
+(R6*P2-R2-^P6)’'-(D3*B2+D3*B1+D3-B3*D2-B3*D1-B3) 
+(R1*P6-R6*P1)*(D3*C2+D3*C1+D3-C3*D2-C3*D1-C3) 

+( R5*P3-R3*P5 ) *( C3*B2-B 3*C2 -C2*B 1 -C2+B2*C 1+B2 ) 

+( R2*P5 -R5*P2 )*( D3’'-B2-B3*D2-D2*B1 -D2+B2*D1+B2 ) 

+( R5*P l-Rl* P5 )*( D3*C2 -C3*D2 -D2*C 1 -D2+C2*D 1+C2 ) 

+( R4*P3 -R3*P4 ) *( C3*B 1 -B3*C 1+C2*B 1 - B2-''C 1 -C 1+B 1 ) 

+( R2*P4 - R4*P2 ) *( D 3*B 1 - B 3*D 1+D2*B 1 - B2*D 1 - D 1+B 1 ) 

+( R4*P1 -R1*P4 ) *( D3*C 1 -C3*D1+D2*C 1 -C2*D1 -Dl+C 1 ) 
+(R3*P2-R2*P3)*(B3+B2+B1+1) 

+(R1*P3-R3*P1)*(C3+C2+C1+1) 

+(R2*P1-R1*P2)-KU3+D2+D1+1) 

=( R6*( P5*D3-'-C2 -P5*C3'>D2+P4*D3*C 1 -P4*C3*D 1 -P3*C3+P2*D3 ) 

+ R5*(-P6-'-D3*C2+P6*C3*D2+P4*D2-'Cl-P4''C2*Dl-P3'-C2+P2*D2) 
+ R4*( - P6'-D3 *C 1+P6*C3"D 1 - P5*D2*C 1+P5*C2*D 1 - P3*C 1+P2*D 1 ) 
+ R3*( P6-'-C3+P5*C2+P4*Cl+P2) 

+ R2*(-P6*D3-P5*D2-P4*D1-P3) ) /( DEN1+DEN2) 

A2=( R6*( -P5*D3*B2+P5*B3*D2-P4"D3*B1+P4*B3*D1+P3*B3-P1--D3) 
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& + R5*( P6-->D3*B2-P6*B3*D2-P4*D2*B1+P4*B2*D1+P3'-B2-P1*D2) 

& + R4*( P6*D3*B1-P6*B3*D1+P5*D2*B1-P5*B2*D1+P3*B1-P1*D1) 

& + R3-’'(-P6*B3-P5*B2-P4*Bl-Pl) 

5 + Rl^'^( P6*D3+P5*D2+P4*D1+P3) )/(DENl+DEN2) 

A3=( R6*( P5^''C3*B2-P5*B3---C2+P4-VC3*B1-P4*B3*C1-P2''-B3+P1*C3) 

6 + R5*( -P6*C3*B2+P6*B3*C2+P4->C2*B1-P4*B2*C1-P2*B2+P1*C2) 

Sc + R4*( -P6*C3*Bl+P6-'-B3-’-Cl-P5*C2*Bl+P5*B2*Cl-P2*Bl+Pr'’-Cl) 

S: + R2*( P6*B3+P5*B2+P4*B1+P1) 

Sc + R1*(-P6*C3-P5*C2-P4*C1-P2) )/(DENl+DEN2) 

IMPEDANCES 



IMP=2*( ( R4*RR+V1*QQ*XK4-V1*PP*XK4**3 )*f D 1*A3+C 1*A2+B 1* A1 ) 
Sc + (R5*RR+Vl-W"XK5-Vl*PP'-XK5**3)-’-(D2*A3+C2*A2+B2’>Al) 

Sc + ( R6*RR+V1*QQ*XK6 - Vl*PP*XK6--'*3 )*( D3->A3+C3’> A2+B3’> A1 ) 

Sc + RR*(R3*A3+R2*A2+R1*A1) 

Sc + V1*Q0*(XK3*A3+XK2*A2+XK1*A1) 

Sc -V1*PP*( K3**3*A3+K2**3-->A2+K1*-’^3*A1 ) ) / ( V1*W) 
IMPR=REAL(IMP) 

IMPI=AIMAG(IMP) 

WRITE( 15 , 1000)F, IMPR, IMPI 

F=F-20 

T1(1)=T1(J) 

T2(1)=T2(J) 

T3(1)=T3(J) 

Q1(1)=Q1(J) 

Q2(1)=Q2(J) 

Q3(1)=Q3(J) 

100 CONTINUE 

1000 F0RMAT(1X,D9. 3,2X,D12. 5,2X,D12. 5) 

CL0SE(UNIT=15) 

STOP 

END 
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APPENDIX B. 
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Figure 49. Real part of the 20" viscoelastic beam waveguide absorber impedances 
at the center of the beam. 
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Figure 50. Iniaginnry part of the 20" viscoelastic beam naveguicle aI)Sorl)er impe- 
dances at tile center of the beam. 
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Figure 51. The driving point frequency response of llie test plate nith a 16" viscoe- 
lastic beam waveguide absorber (dashed) at location 1 and without 
(solid). 
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strained layer beam waveguide absorber (dashed) at location 1 and 
without (solid). 
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Figure 54. Modal damping factors vs. frequency of the test plate \vitli a 20" con- 
strained layer beam waveguide absorber at location 1 and without. 
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(solid). 
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P'igure 56. Modal damping factors vs. frequency of tlie test plate Mitli a 16" viscoe- 
lastic beam waveguide absorber at location 2 and without. 
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without (solid). 
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Figure 58. Modal damping factors vs. frequency of the test plate uitii a 16' con- 
strained layer beam waveguide absorber at location 1 and without. 
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Figure 59. The driving point frequency response of the test plate \>itli a 20" con- 
strained layer beam waveguide absorber (dashed) at location 2 and 
without (solid). 
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Figure 60. Moilal damping factors vs. frequency of the test plate uith a 20" con- 
strained layer beam waveguide absorber at location 2 and without. 
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(solid). 
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MODAL DAMPING VS. FREQUENCY 




Figure 62. Modal damping factors vs. frequency of the test plate uith a 16" viscoe- 
lastic beam waveguide absorber at location 3 and without. 
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without (solid). 
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Figure 66. Modal damping factors vs. frequency of the test plate nith a 20" con- 
strained layer beam waveguide absorber at location 3 and without. 
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